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ABSTRACT 


This volume of the final report on the unmanned Multiple 
Exploratory Probe System (MEPS) details all calculations, deri- 
vations, analyses, and computer programs that support the infor 
mation presented in the first volume. 
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INTRODUCTION 


In Volume One of the report on the Multiple Exploratory 
Probe Syatem (HEPS) the final design Is presented. However, In 
most cases the reasoning or rationale behind many of the design 
decisions are not given to the reader. Volume Tvo alleviates 
this problem by presenting calculations, derivations, computer 
programs, and additional arguments for the final design. 

Several areas are discussed in this volume. First, the 
structural mass calculation and the structural analysis are 
presented. The calculation of the propulsive burns for the 
Earth-Mars transfer are shovn, as veil as the burns required for 
the aerobraking and the satellite. The secondary and main pro- 
pulsion systems are studied to obtain the mass of propellant 
(oxidizer and fuel) required for the entire trip and the size of 
the propellant tanks; engine analysis is also presented in more 
detail in this volume. The necessary equations for the aero- 
braking program are derived, and the lander system is analyzed 
for determination of mass and stagnation temperature. In 
addition the recovery system is optimized. Appendices present 
various programs and supplemental plots. 
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STRUCTURAL HASS CALCULATION 


The total structural mass was calculated by summing the mass 
of each component of the structural system. The structural mass 
will be divided into stringers, bulkheads, the cylindrical 
shells, the caps on the ends of the modules, and the connectors. 
Aluminum will be used for the structural material (specific 
weight - 173.4 lb/ft* ). 

Stringers 

There are 36 longitudinal stringers arranged circumfer- 
entially along the length of each module. Each stringer is 
assumed to have a cross-sectional area of 1 in* * 0. 006944 ft* . 
Therefore, the total stringer mass is 

36 x (module length) x (0.006944 ft*) x (173.4 lbm/ft* ) 
mass (lbs) * 43.365 x length 
Bulkheads 

The bulkheads are I-beams located along the interior circum- 
ference. The cross-sectional area of this beam in 1.25 in* * 

0. 008681 ft* . For the mass of each bulkhead, 

(2* n ) x (12.5 ft) x (0.008681 ft*) x (173.4 lbm/ft» ) 
mass » 118.261 lbm (per bulkhead) 

The (2*n*12. 5 ft) term is the circumference of the bulkhead, 
where 12. 5 ft is the radius of the module. 
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CvAindrAffil ShttlAa 


The structural mass of the cylinder obviously depends on the 
length of the nodule. Using a thickness of 1/4 In (.020833 ft), 
(2*11*12. S ft) x (0.020833 ft) x (173.4 lbm/ft* ) x length 
aass (lbs) * 283.83 x length 
End Caps on Cylinders 

The end cap Is a plate aodelled on the end of the module, 
and two caps are designed for each module. Using the same thick- 
ness as the cylindrical shells (.020833 ft), 

2(it* (12. 9 ft)* x (0.020833 ft) x (173.4 lbm/f t» ) 3 
mass « 3547.98 lbm 
Connectors 

Each connector Is one Inch thick and are five feet In length 
(with the exception of the aerobrake connector, which Is 10 feet 
long). The design of the HEPS vehicle requires a total of four 
5-foot long and one 10-foot long connectors; only two 5-foot 
long connectors and the 10-foot long connector will remain on the 
vehicle during aerobraklng. For the total mass of the connec- 
tors, 

(2* n * 12. 5 ft) x (0.08333 ft) x (length) x (173.4 lbm/ft* ) 
mass for Earth configuration (lbm) * 1135 x 30 ft * 34050 
mass for Mars configuration (lbm) * 1135 x 20 ft * 22700 
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Total Vehicle Hass 


Using the individual calculations given above, the total 
aass of the HEPS vehicle can be deterained: 


Equatorial Lander: 


stringers 

(43.37 lba/ft • 30 ft) 

m 

1305 

lba 

bulkhead 

(4 • 118.26 lba) 

w 

475 

lba 

cylinder 

(283.83 lba/ft • 30 ft) 

a 

8515 

lba 

caps 


■ 

3550 

lba 

total sass 


* 

13845 

lba 

Satelllte/CIC: 

stringers 

(43.37 lba/ft " 25 ft) 

a 

1085 

Vba 

bulkhead 

(3 • 118.26 lba) 

8 

355 

lba 

cylinder 

(283.83 lba/ft • 25 ft) 

s 

7100 

lba 

caps 


8 

3550 

lba 

total aass 

Secondary Propulsion: 

8 

12100 

lba 

stringers 

(43.37 lba/ft • 10 ft) 

s 

435 

lba 

bulkhead 

(2 • 118.26 lba) 

8 

240 

lba 

cylinder 

(283.83 lba/ft ' 10 ft) 

S 

2840 

lba 

caps 


8 

3550 

lba 

total aass 


8 

7065 

lba 

Main Propulsion: 

stringers 

(43.37 lba/ft • 70 ft) 

8 

3040 

lbm 

bulkhead 

(5 • 118.26 lba) 

8 

595 

lbm 

cylinder 

(283.83 lba/ft • 70 ft) 

8 

19870 

lbm 

caps 


8 

3550 

lba 

total aass 



27055 

lba 

Polar Lander: 

stringers 

(43.37 lba/ft • 50 ft) 

8 

2170 

lbm 

bulkhead 

(6 • 118.26 lba) 

8 

710 

lbm 

cylinder 

(283.83 lba/ft ‘ SO ft) 

8 

14195 

lbm 

caps 


8 

3550 

lbm 

total aass 

Total Vehicle (Structural Hass): 

8 

20625 

lbm 

Earth departure « 

111, 500 

lbm 

Hars arrival 

8 

67, 100 

lbm 
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PROPOSED STRENGTH ANALYSIS 


The following analysis approach will bs used to find (a) 
the proper saterial for the stringers, nodule skin, and bulk- 
heads; (b) the correct saterial and thickness of the connectors; 
and (c) the saterial cosposition and overall nusber of pins for 
nodule connections. 

A cosplete finite elesent sodel for the HEPS vehicle is 
being placed on HSC/pal. This sodel includes the proper lengths 
of the nodules and connectors (using alusinus as the initial sa- 
terial for all cosponents) and the payload sass inside each 
nodule. The force cases for dynanic analysis of the systea are 
being obtained fron the propulsion and orbital insertion 
analyses. These cases will include thrust fron the initial de- 
parture burn, acceleration and thrust fron the deceleration burn 
to insert into Martian orbit, and the acceleration and drag from 
the naxinun-force aerobraklng pass. The sodel and forces will 
be translated into a NASTRAN input file using a program available 
on MSC/pal2. 

NASTRAN will be executed using the input cases outlined 
above and the output will be evaluated; to utilize this evalu- 
ation, the following procedure should be applied. If the evalu- 
ation shows that the structural Integrity of the connector is in 
doubt, two cases should be run. First, use aluminum for the 
nodule and a netal natrix composite (HHC) for the connector; 


S 


second, us* ilu*inu« for both components but increase the 
thickness of the connector. If the evaluation shows question in 
the module's strength, the following runs will be conducted se- 
quentially until a proper solution is reached : < 1 ) increase the 

stiffener and bulkhead sizes (areas); (2) change the material of 
the module to MMC and use the original stiffener and bulkhead 
sizes alone; (3) use aluminum for all components but increase the 
thickness of the module; and (4) increase the sizes of all com- 
ponents while still using aluminum. 

Once the proper sizes for a minimum stress on the entire 
vehicle have been determined, an analysis technique outlined by 
R. E. Peterson in Stress Concentration Design Factors (Wiley 
Press, 1974) will be used to determine the stress in a pin hole 
of a connector. Prom this information, strength or failure of 
the connection can be determined. If the pin shows failure, then 
the material of the pins must be changed; however, if the pin 
hole indicates failure, then the thickness around the hole will 
be increased. Also from this stress information, and some 
assistance from Or. W. A. Poster, Jr. of Auburn University, the 
minimum amount of pin connections can be determined. 
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ANALYSIS OF EARTH-MANS TRAJECTORY 


After MEPS has been eoved to the ecliptic plane, a 
propulsive burn will be conducted to start the vehicle on the 
journey to Mars. A Hohsann (sinisus energy) transfer vill be es 
ployed to save fuel, although the tise of flight vill be ex- 
tended. This section vill introduce the analysis of the Earth- 
Mars trajectory, including the sagnitudes of the required propul 
sive burns and the determination of the tine of flight. 

Earth Departure 

The seni-najor axis of the transfer is defined as 

a T « ♦ rg> » *«<4. 908 Ell ♦ 7.477 Eil) ft 

a, « 6.193 Eli ft 

To determine the propulsive burns for the start of the 
transfer, the approach to the problem must be considered. The 
transfer betveen Earth and Mars vill require determination of 
relative velocities, as the situation is not a simple transfer 
betveen tvo orbits about the same body. Nov tvo bodies must be 
taken into account--the Earth and the sun. 

The velocity of the vehicle relative to the sun can be ex- 
pressed in terms of velocities about earth: 

- v © ♦ v -.'© - y 9 ♦ 

vhere V» is the hyperbolic excess speed. This speed can be ex- 



A fundamental aquation used in aatrodynaaics ia tha Via- Viva 
aquation. Thia aquation allova tha calculation of a velocity at 
a point in an orbit if tha paraaatara of tha orbit are knovn: 

V • 

vhara p ia tha gravitational parameter of tha body (planet) vhich 
influancaa tha vehicle, r la tha diatanca from tha body where the 
propulaive burn ia applied, and a ia the semi-major axis of the 
tranafer ellipae. 

Applying the Via-Viva equation to find the velocity of the 
vehicle relative to the aun, the required inputa are 
p 0 « 4. 687 E21 ft* /aec* 
r^ « 4.908 Eli ft 
a - 6.193 Ell ft 
The reaulting velocity ia 

V. # / 0 » 107,383.46 ft/sec 

The velocity of the Earth ia calculated by aaauming the 
Earth ia in a circular orbit about the aun. Uaing the equation 
for velocity in a circular orbit. 



and the appropriate valuea given above, 
v © * 97,722.64 ft/oec 

The hyperbolic exceaa speed can now be determined: 
* V, • /g V^ 

» (107,383.46 - 97,722.64) ft/sec 
« 9660.82 ft/sec 
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The burnout velocity is expressed ss 


V.. - Jv* - |fi» 

At sn radius of 22,567,193.75 ft (3714.153 n si) from the center 
of the Earth, vith 

* 1. 408 E16 ft* /sec* 
the burnout velocity has the value of 
> 36621.86 ft/sec 

The velocity of the vehicle relative to the Earth is given 
by the expression for circular velocity, using the gravitational 
parameter for the Earth and the radius of 22,567,193.75 ft: 

V.., w ■ 24,978.28 ft/sec 

With this value and the value for the burnout velocity the re- 
quired propulsive burn can be calculated: 

AV » V*. - V..,* 

AV « 11643. 58 ft/sec 

The equations for analyzing the propulsive burn to allow 
capture by Mars are similar to those used for the Earth departure 
analysis. 

The hyperbolic excess speed is expressed as 

l v «l * l v --'© * v 6 1 

The velocity of the spacecraft relative to the sun is determined 
by the Vis-Viva equation, vith the following inputs: 
p 0 > 4.687 E21 ft* /sec* 
r£ * 7.477 Ell ft 
a T « 6. 193 Ell ft 
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Substitution into the Vis-Viva squat ion yislds ths valus of this 
vsloclty : 

v .., 0 ■ 70,489.39 ft/sac 

Ths vsloclty of Mars about ths sun is calculatnd under ths as- 
susption that Mars is moving in a circular orbit about ths sun: 



V£ « 79174.22 ft/ssc 

Thus ths hypsrbolic sxesss spssd has ths valus of 
V. ■ 8684.83 ft/ssc 

Ths magnltuda of ths propulsivs burn rsqulrsd for Hars cap- 
turs is sxprssssd as 

&V * V« • /£ — V, , , 

vhsrs V, (r is ths vsloclty of ths spacscraft at ths point of 
clossst approach to Hars. Bscauss ths vshicls is on a hyperbolic 
approach to Hars ths vsloclty of ths vshicls rslativs to Hars is 
sxprssssd vlth ths sama aquation as ths burnout velocity used for 
the Earth departure : 



p£ * 1.S066 E13 ft* /sec* 

r ? * 12,774,573.5 ft (2102.46 n mi) 

V..,^ * 17,643.74 ft/sec 

The elliptic orbit of the vehicle after Hars capture is de- 
fined as 1,640,500 ft (270 n ml) x 108,151,603 ft (17800 n mi). 
Given the radius of Hars, 

rg = 11,134,073.5 ft 
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the aml-Mjor axis of the ellipse can be calculated using the 
f olloving : 

a ■ K-C( 1,640, 300 ft ♦ r > ♦ (108,151,603 f t ♦ r > 3 
a - 66,030,125 ft (10,867.37 n si) 

At periapsis the velocity of the vehicle is deterained 
(using the Vis-Viva equation) to have the following value: 

V,., * 14,596.52 ft/sec 

Using the expression for the propulsive burn required for Mars 
capture, 

AV « V..,g - V,., 

AV ■ 3047.22 ft/sec 
Time of Flight 


The tine required for the transfer fron Earth to Mars can be 
calculated fron the period of the elliptical orbit. The orbit 
period is defined as 


T * 


2n 


»a t • 


Substitution of the appropriate values yields 


p 0 « 4.687 E21 f t* /sec* 
a T - 6. 193 Ell ft 


*r ■ 44, 728, 466. 79 seconds * 517. 7 days 
Since this period is the tine of flight for an entire elliptic 
orbit, the required tine to reach Mars is one-half the period, or 
*T« m a 258. 85 days 
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PROPELLANT ANALYSIS FOR THE SECONDARY PROPULSION SYSTEM 


The Mcondary propulsion system will bs employed upon 
approach to Hars. The purpose of this engine system is to slow 
down HEPS to obtain an elliptic orbit about Hars and to help in 
the final stages of orbit circularization. The analysis pre- 
sented in this section concerns the calculation of the propellant 
mass (oxidizer and fuel) for each AV burn and the required volume 
of the fuel tanks. 

The masses of each HEPS module which will be placed into 
orbit about Hars are presented below. These values do not 
include propellant. 

Aerobrake 12,000 lbm 

Equatorial Lander/Rover 65, 000 lbm 

Satellite System 3,500 lbm 

CIC 3, OOO lbm 

Structure 67, 100 lbm 

The total mass of HEPS excluding propellant is 140, 000 lbm. 

Four AV burns will be required during the circularization 
process at Hars. The first burn will place the HEPS vehicle 
system into a highly elliptic orbit about Hars. During the 
appropriate orbit a second AV burn will be performed at the orbit 
apoapsis to lower the perlapsis into the Hartian atmosphere. 
Following aerobraking, a third burn moves the perlapsis out of 
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the atiotphtr*, and the fourth &V burn vill provide final orbit 


circularization by adjusting tha apoapsia. Tha calculatad values 
of thaaa burns vill nov ba presented. 


Burn 1 (orbit capture) 
Burn 2 (lover periapsis) 
Burn 3 (raise periapsis) 
Burn 4 (adjust apoapsis) 


3047. 90 ft/sec 
75. 1214 ft/sec 
301.6568 ft/sec 
12. 1129 ft/sec 


The secondary propulsion system vill use three engines 
similar to the Space Shuttle Orbiting Maneuvering System (OHS). 
These engines have a specific thrust of 280 seconds. With the 
help of an equation relating the burn to the specific thrust and 
Initial and final masses, the propellant mass (initial mass 
before burn) required for each burn can be calculated. 

&V * Isp*g a *ln(H t /Hr ) 

Mi - H, ..“V'lT-fc > 

Use of this equation vill begin the propellant analysis 
required for the secondary propulsion system; for each burn the 
mass of the propellant (oxidizer and fuel) must be determined. 
The final vehicle mass is 140, 000 lbm. Substitution of this 
mass and the value for the &V burn (12.1129 ft/sec) results in 
the total vehicle mass prior to apoapsis adjustment (or, fol- 
loving periapsis raising): 


M.4 * Hr > 


( 140, 000 lb- > • . < 12 ’ 1 ^ < 280 • 32. i74 ) ) 


Hr a « 140, 188. 36 lbm 
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The propellant naan required for apoapals adjustment la deter- 
mined simply by subtracting the total mass following adjustment 
from the mass prior to the maneuver: 

H,4 « (140,188.36 - 140,000) lbm 
- 188.36 lbm 

Following this procedure, the propellant mass breakdown is given 
in the accompanying table. 


Propellant Hass Required For AV Burns 


Burn 

Humber 


&V Total Vehicle Propellant Hass 

< ft/sec) Hass (lbm) < lbm) 


1 

2 

3 

4 


3047. 90 
73. 12 
301.66 
12 . 11 


203, 026. 33 
146, 173. 91 
144, 962. 03 
140, 188. 36 


58, 8S0. 44 
1213. 86 
4773. 69 
188. 36 


Note that the total vehicle mass on approach to Mars is deter- 
mined to be 203, 026. 33 lbm. The total mass of the propellant 
used during the AV burns is 63, 026. 33 lbm. 

To calculate the mass of oxidizer (N« 04 ) and fuel (HHH) the 
oxidizer/fuel ratio (1.63) will be used. Every 1.63 parts of 
oxidizer is accompanied by 1 part of fuel, for a total of 2. 63 
parts of propellant. From the ratio, 

mass of N«04 * ( 1. 63/2. 63) *Hp 
mass of HHH > < 1. 0/2. 65) *Hp 
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For the total aaaa of propellant given above <63,030 lba) , 
ease of N» 0* * 40, 491 lba 
aass of HHH - 24,540 lba 

The voluaea of the nitrogen tetroxide and aono-aethyl hydra- 
zine are calculated vith the equation for density: 

Voluae * aass/density 

where 

* 85.50 lba/ft* 
f M M H - 33. S3 lba/ft* 

Using the oxidizer and fuel aasses given above, the respective 
voluaes are 

■V hli04 * 473.57 ft* 

* 455.88 ft* 

The shape of the tanks can now be deterained. If cylin- 
drical tanks <25 ft. disaster) are used, the length say be cal- 
culated using 

JL ■ V/<n*R» ) 

For the nitrogen tetroxide, 

“ 0.965 foot 

and for the aono-aethyl hydrazine, 

A Mmh * 0. 929 foot 

Note that the required length is only one foot, which is very 
iapractical. 

Spherical tanks will now be considered. For the radius of 
the tank. 



Calculation of tha aphara radii lor tha oxldizar and fuel tanka 
yields, respectively, 

feet 

R Mmh * 4.774 faat 

Froa thla analyala tha apharlcal tank la tha optlaua daalgn. Tha 
tanka can ba contalnad alda-by-alda within tha cylindrical coa- 
partaant of tha HEPS vehicle; tha radii of tha apharaa may be ln- 
craaaad to five faat for aaaa of construction, thua providing for 
a coapartaant 10 faat In length with sufficient rooa for tha 3 
anginas. 
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ANALYSIS OF THE RAIN PROPULSION SYSTEM 


The following section contains the analysis of the engines 
considered for HEPS. Four engines were cospared on the basis of 
thrust * specific thrust* weight* and burntime; the Space 
Transportation Main Engine (STHE) was chosen for the sain propul- 
sion systes. The necessary data (propellant sass and volume* 
sodule length/tank size) is presented for the STME* and a staging 
analysis is shown. 

Engine Cosoarisons 

Four engines were cospared in the analysis of the sain 
propulsion systes — the J-2* RL-iO-A-i, Space Shuttle Main Engine 
(SSHE)* and the Space Transportation Main Engine (STME). The J-2 
was used for the third stage of the Saturn rockets. The original 
RL-10 engine was used for the early Saturn rockets* and has seen 
use on the Titan; the RL-10-A-1 is sore of an engine design as 
this engine has not been produced. The SSHE is currently in 
operation on the Space Shuttle orbiters* while the STHE is a 
second generation SSME-based engine which also has not gone into 
production. 

By using Newton's Second Law* and assusing the initial mass 
(engine and propellant) to be equal* the four engine candidates 
can be compared: 

£ F » T ■ m«a 

T - ^ 
dt 
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Rewriting the latter equation as an expression lor time. 


dt ■ 
At » 


H»dV 

T 

H«AV 

T 


This final equation is used for calculation of the burn time of 
each engine. Note that the comparisons were made on the basis 
of thrust levels of approximately equal magnitudes; for approxi- 
mately 450, 000 lb of thrust the appropriate number of engines 
must be considered. 

J-2 Thrust 200, 000 lbf x 2 engines 

Isp 416 seconds 

Hass 3460 <6960) lbm 

Burn time 635. S3 seconds 


<702566.2 lbm) • <11641. 26 ft /sec) 

* “ 2 7 <206500lbf ) • <32ri74"Ibm-ft/Ibf-sec* ) 

total mass leaving orbit * 702566. 2 lbm 
delta-V burn required ■ 11641.26 ft/sec 


Thrust 

Isp 

Hass 

Burn time 

SSME Thrust 

Isp 

Hass 

Burn time 

STHE Thrust 

Isp 

Hass 

Burn time 


15. 000 lbf x 29 engines 
433 seconds 

296 <6642) lbm 
564. 40 seconds 

470. 000 lbf x 1 engine 
433 seconds 

6700 lbm 
540. 90 seconds 

435. 000 lbf x 1 engine 
449 seconds 

7455 lbm 
564. 40 seconds 


A short duration burn time is desirable because of the 


decreased risk of course deviation during the burn < ref 3 ) . The 
shortest burn time is achieved by the engine with the greatest 
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thrust) by the above data this engine is ths SSHE. The next 
lowest burn tiss occurs with ths RL-10-A-1 snginss and ths STHE. 
Ths RL-10-A-i was found to bs unfsasibls sines 29 units ars 
rsquirsd to obtain a cosparabls thrust lsvsl. Although ths STME 
has a longsr burn tiss than ths SSHE, ths STHE is dssignsd to be 
sore reliable and less expensive than ths SSHE (ref 2); thus the 
Space Transportation Hain Engine is selected over the Space 
Shuttle Hain Engine. 

Cosparison of the STHE to the J-2 engine is based on thrust, 
weight, burn ties, and design. The two J-2 engines produce 
400,000 lbf of thrust and weigh 6960 lbs; the STHE weighs 
slightly more but produces greater thrust ( 435, 000 lbf ) . The 
burn time of the STHE is considerably less than that of the J-2. 
In addition, the STHE is being designed specifically for reusabi- 
lity and space applications (one design of the STHE nozzle ex- 
pands the flow at the exit to the optimum pressure for operation 
in the vacuum of space). 

STHE Engine Information 

The engine data required for the analysis of the HEPS 

mission will now be presented. Some of the engine particulars 

have been previously stated. 

Isp * 449 seconds 
Thrust » 435,000 lbf 
Hass > 7455 lbm 
Oxidizer/Fuel Ratio * 6. 0 
Area Ratio * 55/141 


19 


The initial of tha HEPS vehicle prior to leaving Earth 

aay ba calculatad uaing tha following analysis: 

AV » Isp*g« *ln(Hi / H, ) 

M. - *... <AV/< I-P-®->> 

Tha vahicla sass bafora tha anginas and tha polar landar systaa 
ara ralaasad is 361,305 pounds (saa pravious aaction). For a 
raquirad delta-V burn of 11641.26 ft/sac to bagin tha Earth-Hars 
transfer. 

Mi - 081,505 lb) ..<«M 1 .2S/<M4 9 , .32.1 7 .». 

Hi • 854,020 lb 

Tha sass of tha propallant is tha diffaranca batwaan tha initial 
mass (vahicla plus propallant) and tha final sass (vahicla only). 
For tha givan conditions, 

H, 1 Hi - H, 

M, « 654,019.4 lb - 361,501.4 lb 
H, « 472,518 lb 

Tha oxldizer/fuel ratio for tha angina is givan as 6. O. 

Evary six parts of liquid oxygan must ba accompanied by one part 
of liquid hydrogen; thus a total of seven parts of oxidizer and 
fuel will be available. Using this development, the masses of 
the liquid oxygen and liquid hydrogen can be determined. 

Hco* » 6/7 - H, * 6/7 - 472,518 lb * 405,015.43 lb 
Mi H« a 1/7 • M # * 67,502.57 lb 
Using the densities of the oxidizer and fuel, and the rela- 
tionship between density and volume, the volumes of the liquid 
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oxygen and liquid hydrogen can be calculated. Knowing these 
voluees will allow the sizing of the fuel and oxidizer tanks. 

• 71.07 lb/ft a at -297° F 
fi.i » 4.42 lb/ft* at -423° F (5*4-23) 

The density is defined as the ease per unit voluee. Therefore. 
•Vi... » 405.015 lb • (1 ft* /71. 07 lb) - 5699 ft* 

* 42.727.7 gallons 

*¥i.„. « 67.503 lb • <1 ft*/4. 42 lb) » 15.271.95 ft* 

*¥iw« » 114.499.6 gallons 

A pressure vessel- is noraally spherical, or cylindrical with 
heslspherical ends. The diaseter of the HEPS vehicle must be 
considered to detereine which type tank will hold the liquid 
oxygen and hydrogen. For the diaaeter of 25 feet, the voluae of 
a spherical tank is 

^*4/3- w <12. 5 ft)* ■ 8161.23 ft* 

This voluae falls between the required voluaes for the oxidizer 
and fuel. Thus, a spherical tank will be eaployed for the liquid 
oxygen, and the cylindrical /heaispherical tank will be used for 
the liquid hydrogen. 

For the calculated voluae of liquid oxygen the corresponding 
tank size is deterained to be 

voluae « 5699 ft* * 4/3* n*R* 

R » 11.08 ft 

If boil-off of the liquid oxygen (12.0 lb/hr) is considered, the 
actual size of the tank aust be increased to account for the 
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expendable oxidizer. Based on a thirty day transport and con- 
struction period, the voluae of the L0« lost to boil -off is 
aass > (720 hrs>‘<12.0 lb/hr > « 8640 lb 
voluae > (8640 lb) Ml ft*/71.07 lb> - 121.57 ft* 

This additional voluae yields an increase in the diameter of the 
spherical tank to 22. 5 feet. 

The voluae of liquid hydrogen is auch larger than that of 
the liquid oxygen and, as mentioned previously, a cylindrical 
tank vith hemispherical endcaps will be required. For the tank 
to fit snugly Inside the HEPS vehicle (diameter of 25 feet), the 
length of the tank can be calculated: 

voluae « 4/3- *•( 12. 5)* «■ nM12. 5)»*h » 15,271.95 ft* 
length ■ h ■ 14. 4 ft 

The tank size vill increase under consideration of boil-off. The 
rate for LH« is 18. 0 lb/hr, and for the same thirty day period 
used earlier, 

aass ■ (720 hrs)M18. 0 lb/hr ) * 12,960 lb 
voluae > (12,960 lb) Ml f t* /4. 42 lb) * 2932.13 ft* 

The change in the length of the tank is nov determined: 

length « h ■ (15,271.95 ♦ 2932.12 - 8181. 23) /490. 87 
h » 20.45 ft 

Since the endcaps have a radius of 12. 5 feet, the total length of 
the LH. tank is 

25 ft ♦ 20.45 ft * 45.45 ft 
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If the two tanks art mounted bulkhead to bulkhead, the total 


length of the main propulsion module Is 
45.45 ft ♦ 22.5 ft « 67.95 ft 
Staging Analysis 

In order to determine the optimum number of stages for Earth 
departure, the fallowing equation is used: 

x- (y - s] n 

The results of the analysis performed on the STME are plotted 
below. Prom this plot, one stage is shown to be the optimum con- 
figuration for the engines. 


PLOT TO OPTIMIZE STAGING 
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AEROBRAKE ANALYSIS 


The atrobrake will be uaed for orbit circularization about 
Mara. Thia aaction will praaant tha dataralnatlon of tha weight 
of tha aarobraka and tha calculation of tha propulalva burna re- 
quired for tha circularization analyaia. 

Wflah* fitlREBinRiian 

Uaing inforaation provided by Bill Wlllcockaon, OTV Program 
Manager at Martin Marietta Aaroapaca * Denver ), tha maaa of an 
aarobraka can be alzad with tha aarobraka area. From valuaa pre- 
aantad in reference 7, 

Area - 142 ft* 

maaa of rigid aurfaca inaulation (RSI) » 401 lbm 
maaa of flexible aurfaca inaulation (FSI) * 8890 lbm 
atructura weight * 11032 lbm 

Tha weight of tha RSI will remain 401 lbm ainca a diameter of 25 
feat ia uaad for tha Martin Marietta brake aa wall aa tha 
propoaad brake. Tha weight of tha FSI will require a calcula- 
tion. For tha Martin aarobraka, 

A * ; C < 142 ft)* - (25 ft)*] * 15345.895 ft* 

4 


Obtaining a weight to area ratio, 

weight 8890 lbm 

area * 15345.895 ft* 

For tha HEPS aarobraka. 


.57931 Ibm/ft* 


A « ; c <93 ft)* - (25 ft)*) = 6597.345 ft* 

4 
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The weight of the FSI is calculated using the vsight/area ratio 
previously determined: 

W# • i - < . 37931 lbs/ft* ) • < 6397. 343 ft* > 

- 3821.893 lba 

To detersine the weight of the aerobrake structure, size the 

veights using a weight /area ratio: 

» . 11032 lba 

^ • <93 ft)* %• <142 ft)* 

* 4937.701 lba 

Propulsive Burns Used For Mars Orbit Circularization 

Although aerobraking will be applied during the HEPS 
Mission, couplets orbit circularization will require propulsive 
burns using the secondary propulsion systen. These burns must be 
considered for three different Maneuvers: lowering the periapsis 

prior to aerobraking, and raising the periapsis and adjusting the 
apoapsis after braking. Calculation of each V will be made by 
applying the Vis- Viva equation: 

V » * <2/r - 1/a) 

where r is the length of either the periapsis or apoapsis, 
measured from the center of Hars, and a is the semi -major axis of 
the elliptic orbit. 

To decrease the periapsis, the burn will be applied at the 
apoapsis of the initial elliptic orbit about Hars. The desired 
periapsis altitude has been determined to be 314,976 ft <51.84 
nautical miles). The lengths of the initial periapsis and apo- 
apsis are 120,398,076.5 ft <19848.27 n mi) and 12,774,373.3 ft 
< 2102. 46 n mi ) , respectively, measured from the center of Hars. 
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The lengths of the MRi“»«jor axes of th* two different 
•lliptical orbits (sans initial apoapsis, two different peri- 
apaaa > ara datarainad t 

r p « 12,774,373.3 it (2102.46 n ai) 

a T| - K(r* ♦ r # ) ■ 66,666,323 ft (10,973.37 n ai) 
r, » 11,449,049.5 ft (1664.31 n ai) 

•r z » 66,023,563 ft (10,666.29 n ai) 

Applying tha Via-Viva aquation, tha valocitiaa at tha apoapsis 


for aaeh alliptical orbit ara oalculatad: 

1 [ e,A=> sac '120596076.5 ft 

V« ■ 1546. 98 ft /sac 


1 

66666325 ft 



Vi ■ f( 1.5066 E15 . < 2 — . — , - )~1 ^ 

* px.ouoo 613 B#c ' 1 120598076. 5 ft 66023563 ft J 

V. > 1471.83 ft/aac 

Tha propulsiva burn raquirad to lover tha pariapsis is datarainad 
by taking tha diffaranca between tha apoapsia valocitiaa givan 
abova j 

A V « V, - V, - -75. 1214 ft/aac 

Tha sinus sign indicataa tha burn will ba appliad in tha direc- 
tion opposita that of tha HEPS vahicla (ratrofira). 

Tha aaaa analysis is parforaad for tha burns to raisa tha 
pariapsis and apoapsia. Tha nacassary inputs and output are pre- 
sented : 
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to raiat the pariapala to 1,640,300 ft (270 n ai)t 


r* « 12,717,736.74 ft (2093.11 n mi) — froa 

r, * 11,449,049.3 ft (1664.31 n al) 

a T * 12,083,391.3 ft (1988.71 n ai) 

r, - 12,774,573.3 ft (2102.46 n al) 

a T » 12,746,133.2 ft (2097.79 n al) 

AV * 301.6568 ft/aac 

to raiaa tha apoapala to 12, 774, 373. 5 ft ( 2102. 46 
tha cantar of Mara: 

r, * 12,774,573.5 ft (2102.46 n al) 


r« 

* 

12, 717, 736. 7 

ft 

(2093. 11 

n 

al ) 

at 

U 

12, 746, 153. 2 

ft 

(2097.79 

n 

al ) 

r. 

M 

12, 774, 573. 5 

ft 

(2102. 46 

n 

al) 

a t 

* 

12, 774, 573. 5 

ft 

(2102. 46 

n 

al ) 


AV « 12. 1129 ft/aac 


prograa 


n al ) froa 
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DERIVATIONS FOR THE AEROBRAKINO PROGRAM 


The program included in Appendix A is used to execute the 
iterations for the aerobraking process. With inputs concerning 
the orbit of a vehicle about Mars, and parameters of the aero- 
brake, the complete aerobraking passage can be analyzed. The 
output presents the time for aerobraking, the drag forces that 
act on the aerobrake, and the parameters of the final orbit. 

The program requires several derivations — the location of 
the intersection of a circle (Mars atmosphere) and an ellipse 
(vehicle orbit); the length of segment between the intersection 
points (total distance travelled within the atmosphere); and the 
drag coefficient of the aerobrake. 

InttrftoUon Point# 

The equations of an ellipse and a circle are given, respec- 
tively, as 

$ * t • * 

(x - ae)* ♦ y* ■ r* 

where a is the semi -major axis and e is the eccentricity of the 
orbit; and a*e is the location of the center of the circle repre- 
senting the Martian atmosphere (i.e. , the center of Mars). In 
addition, the trajectory equation, which gives the location of 
any point on the ellipse, is defined as 

P 


r * 


1 ♦ e • cos V 0 
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In the above aquation all of tha variables (semi-latus rectus, 
eccentricity, and argusent of periapsls) are given parameters of 
an elliptic orbit. 

Solving for y* in the two equations. 


y* ■ 
y* - 

Equating the y* 


b* 


b* x a 

a* 


r a - <x - «•)* 
terms. 


b* 


b a x* 

a* 


r* - (x - ae>* 


b* 



* r * 


x* 


2aex - a*e a 


(1 


b* 

=r> *x» 


2aex ♦ <b a 


r a ♦ a a e a ) * 0 


Using the quadratic equation to solve for x, 

2ae ± v/4a a e a - 4»(b a - r a * a a e a )Ml - b a /a a >' 

2* ( 1 - b a /a* ) 

Choosing only the negative value of the square root (due to the 
geometry of the problem), the x-location of the intersection 
points in known. Substitution of x back into an expression for y 
will yield the complete location of the points. 

Segment Length 

The values of x and y obtained as the intersection points 
will be used in this derivation. The angle created between radii 
from the center of Mars is denoted as 6, s is the segment length, 
R is the radius, c is the chord of the arc, and d is the distance 
from the center of Mars to the x-position of the intersection 
points. 
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Using the geometry of an arc, 


T 
1 

c « 2y 
d * (x - ae) 

R » <d* ♦ y* )H. 

0 » 2 • tan' 1 < jrj* > * 2 'tan 1 (— ^) * 2* tan" 1 ( • ) 

20 a x-ae 

Finally, 

s ■ R0 



Thus the segment is known, and this value can be used to deter- 
mine the change in velocity due to drag forces during aero- 
braking. 


>etermination of the Drag Coefficient 


From hypersonic equations for a cone, the drag coefficient 


is given as 


C» * 2*. sin* 0 V ♦ <1 - 3* sin* 0 v )*sin* a 


where 


0„ » cone half-angle 
a * angle of attack 

For the HEPS mission the design for zero angle of attack (using 
momentum wheels and the cone's inherent stability) allows cancel- 
lation of the second term. Therefore (1:681), 

C» * 2* sin* 0 y 


30 



ANALYSIS FOR THE OBSERVATION SATELLITE 


This section will contain the calculations of ths propulsive 
burns required for orbital transfer of the satellite. In 
addition, the determination of the solar array panel (applicable 
to the CIC as veil) vill also be presented. 

Propulsive Burns 

A transfer between a 1,640,500 ft (270 n si) orbit and a 
2,313,105 ft (327.64 n si) orbit vill be required to put the 
satellite into the observation orbit. For these calculations a 
Hohsann (sinisus energy) transfer vill be assumed. 

Calculation of the altitudes is the first step: 
rg * 11,134,073.5 ft (1632.47 n mi) 
r, « (11,134,073.5 ♦ 1,640,500) ft 
» 12,774,573.5 ft (2102.46 n mi) 
r, ■ (11,134,073.5 ♦ 2,313,105) ft 
> 13,447,176.5 ft (2005.25 n mi) 

The gravitational parameter of Mars is given as 


Mg * 1.5066 E1S ft* /sec* 

For a Hohmann transfer, first calculate the circular 
velocities of the tvo orbits: 


V e 


V, 



10,859.7819 ft/sec 


V, 


1.5066 E15 ft* /sec* 
13447178.5 ft 


lO, 584. 8341 ft/sec 
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Nov determine the seml-ma jor axis of thv transltr ellipse: 
a, « %• <r» ♦ r, ) » 13,110,876 ft (2157.82 n nl) 

To obtain the propulsive burn required to leave the Initial 
orbit, the Vis-Viva equation of astrodynamics (see trajectory 
analysis section) will be used. The result Is 
V T| * 10,998.2401 ft/sec 

The burn is found by subtracting the circular velocity fros the 
velocity at the periapsis: 

^ V. * V T| - V» » 138.4 ft/sec 
The speed at the apoapsls and the propulsive burn required 
to achieve the final circular orbit are calculated in a similar 
manner : 

Vt * 10,448.0164 ft/sec 
& V, - 136.64 ft/sec 

Using these burns the mass of propellant required for the 
transfer can be calculated (the proper equation may be found in 
the section on the secondary propulsion system). First obtain 
the mass ratio for each &V: 

<j| L ) * 1.0155 ( 2“ > - 1.0153 

fir t n, « 

The final mass of the satellite in the observation is approxi- 
mately 3500 lbm. Backing out the mass required by the second AV, 
Ht * H# ♦ M, 

1.0153* Hr - Hr ♦ H # 

0.0153- (3500 lbs) * M, 

H, > 53.55 lbm 
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Obtaining the propellant aaaa uaad In the first AV with the same 


calculation#, 

1* 0155*11, » H, ♦ M # 

0. 0155* H, - M, 

0.0199* <3500 lba ♦ 53.59 lba> * M, 
H, « 55.08 lba 


Tha total aaaa of propellant ia 

M, - M, ♦ H, 

T 

M # * 108.93 lba 
T 

Arti tar th t .SftAiE-Arny 


The area of a aolar array panel ia calculated uaing 


A« 


E$ 

S **| • F • coa T 


where 


(ref 5) 


P* * available power <1500 watta) 

S » aolar intenaity at Mara (54.14 aW/ft* ) (5:410) 


1 


BSX?£-2li£BBi:-2?_§C£S£ o 

power Input of aun * 

aua-total of array deaign and degradation factora 

aiac. aaaeably and degradation 0. 95 

radiation (for ailicon cells) *0.74 

configuration (flat plate array) 1.00 

F » 2.89 (5:123-125) 


r * angle between aun 'a rays and the noraal to the 
panel 

T * 0 

cos r * i 


Calculation of the area yields 


A« 


1500 W 

<54.14 W/ft* )• (0.096) • (2.69) *1 


* 108 ft* 
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CALCULATION OF APPROXIMATE HASS OF MARS LANDER 


To calculate the approximate mama of a Mara landar the 
aaaaaa of tha major componanta of tha landar muat ba aatimatad. 
Thaaa major componanta ara ( 1 > tha Sample Return Vehicle ( SRV > ; 
(2) tha rovers; <3) tha automated laboratory; <4) an upper and 
lover aaroahall; <S) a platform for tha SRV to ait upon; 

<6> landing gear; and (7) a recovery ayatam consisting of a solid 
rocket motor and three parachutes. 

Maas of the Sample Return Vehicle (SRV) 

The mass of the solid rocket booster that will propel the 
SRV can be obtained from a program written to simulate the launch 
of a solid rocket booster; booster specifications include the 
fuel, payload, and planet of launch. This program, titled 
# Stages", can be used to study the effect of changing propellant 
mass on the final altitude and velocity achieved by the rocket. 

To use "Stages” (listed in Appendix A), the following para- 
meters for a launch must be known or assumed: 

1. the combustion temperature of the propellant <°R) 

2. the density of the propellant (lbm/ln 4 ) 

3. the propellant cross-sectional diameter (feet) 

4. the propellant burn rate (in/sec) 

5. the specific heat ratio and perfect gas constant 

for the burning propellant 
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6. the radius of the planet <n. ai. ), tha latitude of 

launch site (degrees), the angular velocity of 
the planet surface, and the gravitational accel- 
eration on the planet surface 

7. the desired altitude after launch 

8. the ease of the final payload to be put into orbit. 

The first seven paraseters were designated for a launch from 

a pole of Mars of a rocket propelled by the propellant DB/AP- 
HMX/AL (Double Base/ Aluminum Perchlorate-Cyclotetramethylene 
Tetranitramine/Aluminum), selected for its high combustion temp- 
erature (6700 degrees Rankins) and burn rate (.95 in/sec). The 
rocket booster vis designed to have a propellant cross-sectional 
area of 2.91667 feet and a deadweight ratio (ratio of booster 
non -propellant mass to total booster mass) of 0. 12. The desired 
orbit vas specified to be circular at an altitude of 270. 0 nau- 
tical miles. 

The eighth parameter (payload mass) vas designed to be a 
lightweight vessel that would carry up to 100 lbm of Martian soil 
and air samples in a refrigerated chamber; on board the ship 
would be a small reaction control system and an aeroshell. The 
mass of this vehicle is estimated to be 1000.0 lbm (200 lbm for 
the refrigeration chamber, 500 lbm for the reaction control sys- 
tem, 100 lbm of samples, 50 lbm for the aeroshell, and 150 lbm 
for onboard guidance and control computers). 

Once these parameters for the SRV launch have been speci- 
fied, a "target burn time” (which is equal to the mass of pro- 
pellant divided by the mass consumption rate) is entered into the 
"Stages" subroutine named "Launch". This subroutine is a numerl- 
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cal integration of tha aquation of aotion for a single-stage 
rocket being launched in a gravity field# and vill fire for the 
entire "target burn tiae" unless the propellant is coapletely 
consuaed or the target altitude is reached. As the launch pro- 
ceeds# the spacecraft is rotated through a "pitch prograa"# 
arbitrarily selected to vary the direction of the rocket's 
weight vector as its altitude increases. 

To optiaize the propellant aass# and thus the initial aass 
of the SRV# five plots of data froa "Stages" are constructed (see 
Appendix B). These plots show: 

1. variation of final altitude with "target burn time” 

(Figure B. 1 ) 

2. variation of final velocity with "target burn time” 

( Figure B. 2 ) 

3. variation of payload ratio (payload aass/initial 

aass) with "target burn tiae” (Figure B. 3) 

4. variation of "excess aass* (aass excluding payload 

aass after launch) with "target burn time" 

( Figure B. 4 ) 

5. variation of final acceleration with "target burn 

time” (Figure B. 5) 

The first plot is used to determine a minimum value for 

« 

"target burn time" (TBT) by observing that below a particular 
value for TBT the desired altitude is not reached due to insuf- 
ficient propellant mass. The second plot is then used to find 
the range of values of TBT above the minimum altitude value for 
which the final velocity is at least sufficient to achieve a cir- 
cular orbit at the design altitude. 
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The third, fourth, and fifth plota ara uaad to find tha op- 
tiaua vaiua of TBT froa tha ranga of TBT valuas datarainad with 
tha firat two plota. An optiaua payload ratio can ba aalactad 
froa tha third plot; an optiaua "excess mass" can ba aalactad 
froa tha fourth plot; and an optiaua final accalaration can ba 
aalactad froa tha fifth plot. 

For tha SRV tha payload ratio vaa optiaizad bacauaa tha 
ratio will yiald a ahortar and laaa aaaaiva booster than tha 
booatar for which tha final accalaration ia ainiaua. Froa Figure 
B. 1 tha daairad altituda ia raachad for TBT graatar than or equal 
to 3S0 aaconda. Froa tha aacond plot tha required orbital velo- 
city ia achieved only for TBTe ranging froa 100 to 460 seconds; 
therefore tha optiaua ranga of TBTe ia between 350 and 460 
seconds. Froa tha third plot tha TBT for tha highest payload 
ratio is found to ba 350 seconds. Froa tha fourth plot axcaaa 
aaaa ia seen to ba a ainiaua for tha optiaua ranga at TBT of 350 
seconds. Froa Figure B. 5 final accalaration is seen to ba a max- 
iaua for tha optiaua ranga at TBT equal to 350 seconds. 

Tha value used for tha payload ratio of tha SRV, baaed upon 
optiaization by use of "Stages*, is 0.066120. This ratio yields 
an initial SRV mass of approximately 14, 700 lbm. 

Haas of the Upper and Lower Aeroshell 

The masses of the upper and lower aeroshell can be estimated 
by determining the approximate geometry of the aeroshells and 
selecting a material with which the aeroshells will be made. The 
material that is selected for tha shall must be strong enough to 
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withstand large aarodynaaic forces and aarotharmodynaaic haating 
incurrad upon daacant through tha Martian ataoaphara. Tha 
aaxiaua taaparatura that will occur on tha landar during daacant 
will ba located at tha stagnation point of tha vassal, which is 
located at tha canter of tha lower aeroahell. Tha stagnation 
point taaparatura that tha landar encounters at an altitude of 
100 nautical alias (calculated as 1630 degrees Ranklna in Appen- 
dix C) is used to dataraina tha type of aatarial to be used for 
tha upper and lower aaroshall (as a first approxiaation ) . 

An "outer blanket" of carbon-carbon heat -transfer resistant 
tiles, or a ona-piaca carbon-carbon sheet, will cover tha bottom 
of tha lower aaroshall. Tha inner part of tha lower aeroshell 
and tha upper aaroshall will ba composed of CLAD 2014 aluminum 
alloy (density of .101 lbm/ in*). Modeling tha upper aeroshell as 
a conic frustrum 23 feet high, with a base diameter of 25 feet, 
a top diameter of seven feet, and a thickness of . 30 inches, an 
approximate upper aeroshell volume of 233, 989. 7 cubic inches and 
a mass of approximately 23, 633. 0 lbm are determined. Modeling 
tha lower aaroshall as a segment of a sphere with a base diameter 
of 25 feet, a segment height of five feet, and a thickness of . 30 
inches yields an approximate lower aeroshield volume of 20, 722 
cubic inches and a mass of approximately 2093. 0 lbm. 

Hass of Platform and Landing Gear 

A metal disk twenty-five feet in diameter and one-half inch 
thickness is used to model the platform which the SRV, rovers, 
and autonomous laboratory sit upon; this platform has a volume 
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of 33, 342 cubic inches. A strong satsrial that can withstand the 
effects of the exhaust pluse of the launching SRV is needed to 
coaprise the platform; AM-355 stainless steel is chosen for its 
favorable resistance to high temperature and corrosion. AM-355 
stainless steel has a density of . 282 lbm/in* , so the mass of the 
platform is approximately 9966 lbm. 

Each strut of the landing gear was modeled as a quarter -inch 
thick AH-3SS stainless steel pipe, one foot in outer diameter and 
five feet long, fastened to a square AM-355 stainless steel pipe 
vith sides four feet in length and a thickness of one-half inch. 
The total volume of each strut is 1857 cubic inches and the 
total mass of each strut is approximately 525 lbm. The landing 
gear system vill consist of four struts so the total landing gear 
mass is 2100 lbm. 

Mass of the Rovers and Automated Laboratory 

Each rovers is to be no more massive than 2500 lbm (5000 lbm 
for the two rovers on each lander). The mass of the automated 
laboratory vill not exceed 1000 lbm. 

Mass of the Lander Recovery System 

The mass of the recovery system for the lander was deter- 
mined by use of a program written by D. Bell. The program cal- 
culates the optimum recovery system mass, consisting of one to 
six parachutes and a solid fuel retrorocket, for a generic ve- 
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hide of ipvcifiibl* mass (Minus rscovsry systsa aass) that is 
landing on the surface of Earth or Kars. The program inputs are 

1. the mass of the vehicle vithout recovery system 

2. the desired terminal velocity for the main para- 

chutes 

3. the number of parachutes desired 

4. the specific impulse, thrust, and mass fraction of 

the solid rocket motor used for descent 

9. the required velocity upon impact with the planet 
surface 

6. the desired height above the ground at which a 

constant -velocity descent of the vehicle begins 
(the rocket is fired such that the thrust equals 
the veight of the vehicle — "constant velocity 
falling height") 

A set of plots can be obtained by making a series of runs of 
this program (see Figures D. 1, D. 2, D. 3). These plots are used 
to optimize the main chute terminal velocity, impact velocity, 
constant velocity falling height, the mass of the parachute sys- 
tem, and the mass of the solid rocket motor required to land the 
vehicle. 

The total (approximate) mass of a lander is found to be 
58, 500 lbm, excluding the mass of the recovery systep. A set of 
runs of the optimizing program were made using this value of the 
vehicle mass; the results can be seen in the figures of Appen- 
dix D. The optimum main chute terminal velocity for this vehicle 
is determined to be 75 feet per second. The optimum impact velo- 
city for the vehicle is 10 feet per second (assuming that the 
terminal velocity under consideration for the vehicle is the op- 
timum value). The optimum constant velocity failing height is 
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five feet (tituiing that the tarainal velocity and iapact velo- 
city are optiaua valuea. The optimum recovery eystem mass for 
the lander (using the aforementioned variables) is determined to 
be as follows: 

parachute system mass * 1610 lbm 

solid rocket motor mass ■ 2639 lbm 

total recovery system mass * 4249 lbm 
A more detailed breakdown of the mass of the recovery system is 
shown in Appendix D. This breakdown is the final output of the 
optimizing program. 


W LI 3 a»x*J -5 J->:f IlI-l m tFXTl WLiLlLm->mA\4m* ix-m 


As stated at the beginning of this section, the approximate 
total mass of the Hars lander is the sum of the masses of its 
major components: 


COMPONENT 

HASS (lbm) 

SRV 

14, 700 

Upper Aeroshell 

23, 633 

Lower Aeroshell 

2, 093 

Platform 

9, 966 

Landing Gear 

2, 100 

Rover Systems 

S, 000 

Laboratory 

1,000 

Parachute System 

1,610 

Solid Rocket Motor 

2,639 

TOTAL APPROXIMATE LANDER MASS 

62, 749 
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APPENDIX A 

Computer Program for the Amrobraking Analysis 
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THIS PROGRAM RUNS THROUGH AN AEROBRAKE ANALYSIS OF MARS 
INPUTS ARE MADE IN SI UNITS, AND OUTPUT IS WRITTEN IN 
ENGLISH UNITS 

REAL MU, HASS, MASSE 

OPEN ( UNIT-7, FILE- ' AEROBRKE. DAT ' , STATUS* ' OLD ' ) 

PRINT*, 'INPUT THE PERIAPSIS ALTITUDE IN km ' 

READ (6,*) PERIAP 

DEFINE THE PERIAPSIS FROM THE CENTER OF MARS 
PERIAP - PERIAP * 3393.5 

PRINT*, 'INPUT THE INITIAL SEMI-MAJOR AXIS IN km ' 

READ (6,*) AINIT 

PRINT*, 'INPUT THE DESIRED APOAPSIS DISTANCE IN km ' 

READ <6, •) APOAPF 

CALCULATE THE TIME FOR THE FIRST 1.5 ORBITS BEFORE MAKING 
PERIAPSIS CHANGE < AT APOAPSIS) FOR AEROBRAKING PROCESS. 

THE ELLIPTIC ORBIT FOR THIS PERIOD IS (500 km X 33363 km) 

SHAJ1 * .5* ((500. *3393. 5) *(33363. *3393.5)) 

THE GRAVITATIONAL PARAMETER MU (kg A 3/sec A 2> 

MU * 42656. 

PI * 3.141592654 

PERD1 * 2*PI*SQRT(SHAJ1**3/MU) 

TIME1 « 1. 5*PERD1 

CALCULATE THE PARAMETERS OF THE ELLIPTIC ORBIT ABOUT MARS 
USING THE PERIAPSIS FROM INPUT. THIS ORBIT IS ACTUALLY ONLY 
HALF AN ORBIT, MAKING THE JOURNEY FROM APOAPSIS TO AEROBRAKING 
PERIAPSIS 
PERIOD * 0.0 
APOAP -0.0 

CALL PARAHS ( PERIAP, APOAP, AINIT, PERIOD ) 

PERDHR • PERI0D/3600. 

DETERMINE THE TIME FOR THE HALF-ORBIT FROM THE APOAPSIS TO 
THE PERIAPSIS OF AEROBRAKING 
TIME2 * . 5*PERI0D 

OBTAIN THE INPUTS FOR THE AEROBRAKING PROCESS 

PRINT*, 'INPUT THE ATMOSPHERIC DENSITY FOR THE ALTITUDE (kg/m A 3> ' 
READ (6,*) RHO 

PRINT*, 'INPUT THE MASS OF THE SPACE VEHICLE (kg) ' 

READ (6,*) MASS 

CONVERT THE MASS TO ENGLISH UNITS 
MASSE * MASS*32. 174/14. 57 

PRINT*, 'INPUT THE HALF-ANGLE OF THE CONICAL AEROBRAKE (deg) ' 

READ (6,*) THETA 

PRINT*, 'INPUT THE DIAMETER OF THE AEROBRAKE (m) ' 

READ (6,*) DIAH 
DIAM * DIAM/1000. 

DETERMINE THE AREA OF THE CONICAL AEROBRAKE 
THETAR * THETA*PI/180. 

PART * 1. /(TAN (THETAR )**2) 

AREA * PI*(DIAM/2. >**2*SQRT(1. *PART) 

DETERMINE THE DRAG COEFFICIENT OF THE AEROBRAKE 
BASED ON NEWTONIAN METHODS 
CD - 2. *SIN( THETAR)* *2 

CALCULATE THE ENGLISH-UNIT COUNTERPARTS OF THE ABOVE VALUES 
PERAPE * PERIAP*3280. 839895 
AINITE - AINIT*3280. 839895 
DIAME - DIAM*3280. 839895 
AREAE * AREA* 10763910. 42 



■ RHOE - RH0*0. 0019435035 

T THE PER I APS IS, APOAPSIS, SEMI -MAJOR AXIS, AND AEROBRAKE DIAMETER 

ARE CONVERTED FROM km TO ft. THE AREA IS CHANGED FROM km* 2 TO ft *2 
AND THE DENSITY FROM kg/m A 3 TO alugm/ft A 3. 

TIMTTL * TIME1+TIME2 
TIMTLH * TIMTTL/3600. 

WRITE (7, *)' AEROBRAKE ANALYSIS ' 

WRITE (7, *) ' ' 

WRITE (7, *) ' ' 

WRITE (7, ») ' ' 

WRITE (7, • ) ' HALF-ANGLE FOR CONICAL AEROBRAKE (deg): ', THETA 
WRITE (7, O' DIAMETER OF THE AEROBRAKE (ft): DIAHE 

WRITE (7, • ) ' SURFACE AREA OF THE AEROBRAKE <ft A 2): ' , AREAE 
WRITE (7, *) ' HASS OF SPACE VEHICLE (lb): MASSE 

WRITE (7, * ) ' ' 

WRITE (7, *) ' ' 

WRITE (7, *>' ATMOSPHERIC CONDITIONS: ' 

WRITE (7, *)' PERIAPSIS FROM CENTER OF MARS (ft): ' , PERAPE 
WRITE (7, *)' DENSITY (slug/ft *3): ' , RHOE 
WRITE (7, * ) ' ' 

WRITE (7, * ) ' INITIAL ORBITAL PARAMETERS: * 

WRITE (7, * ) ' ' 

WRITE (7, 130) 

WRITE (7, • ) ' ' 

WRITE (7, 140) PERAPE, APOAPE, AINITE, PERDHR 
WRITE (7, * ) ' ' 

WRITE (7, * ) ' ' 

WRITE (7, *) ' APPROX TIME (hr) PRIOR TO AEROBRAKING: ', TIMTLH 
WRITE (7, » ) ' ' 

WRITE(7, *)' ' 

WRITE(7, * ) ' ' 

WRITE ( 7, * ) ' AEROBRAKE PROCEDURE : ' 

WRITE(7, *) ' ' 

WRITE (7, 100) 

WRITE (7, 110) 

WRITE (7, * ) * ' 

SET INITIAL CONDITIONS FOR VARIABLES PRIOR TO DO-LOOP 
SMAJ * AINIT 
X * 0.0 
Y * 0.0 
SEG *0.0 
PHI * 0.0 
ASECTR * 0.0 
TIMTTL * 0.0 

DO 50 I * 1,500 

ECCNTY * ( APOAP-PERIAP ) / ( APOAP+PERIAP ) 

SMIN * SORT (SMAJ* *2* ( 1. -ECCNTY**2) ) 

CALCULATE THE INTERSECTION POINTS OF THE ELLIPTIC ORBIT 
: AND THE MARTIAN ATMOSPHERE 

CALL NTRSEC ( SMAJ, SMIN, ECCNTY, X, Y ) 

: CALCULATE THE LENGTH OF SEGMENT OF THE ELLIPTIC ORBIT 

ENCLOSED BY THE MARTIAN ATMOSPHERE 

CALL SEGMNT ( X, Y, SMAJ, ECCNTY, SEG, PHI, ASECTR) 

J DETERMINE THE VELOCITY OF THE SPACECRAFT AT PERIAPSIS 
VELCTY * SORT (MU* (2. /PERIAP-1. /SMAJ) ) 


THE SEMI -MAJOR AXIS IS THE 'OLD* SEMI-MAJOR AXIS 
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CALCULATE THE DRAG ON THE VEHICLE DURING THE AEROBRAKING PROCESS 

UNITS ARE ( KG*KM/SEC A 2 ) AND (LB) 

DRAG « . 5*CD*(RHQ*1. E9)*VELCTY**2*AREA 
DRAGE » . 5*CD»RH0E*VLCTYE**2*AREAE 

DETERMINE THE TIME (IN MINUTES) OF THE AEROBRAKE PASSAGE 
TIME * 2. *ASECTR*SQRT(SMAJ/HU)/SMIN 
TIME « TIME/60. 

DETERMINE THE NEW SEMI -MAJOR AXIS 
ENRGY1 * -MU/(2. *SMAJ) 

SHAJ » -MU/ (-2. *DRAG*SEG/HASS*2. *ENRGY1) 

DETERMINE THE PARAMETERS OF THE NEW ELLIPTIC ORBIT 
CALL PARAHS ( PERI AP, APOAP, SHAJ, PERIOD ) 

PERIOD OF THE ORBIT IS IN HOURS 
PERDHR * PERI0D/3600. 

CONVERT SI UNITS TO ENGLISH UNITS 
SMAJE > SHAJ *3280. 839895 
APOAPE > APOAP *3280. 839895 




CHECK IF THE APOAPSIS IS LESS THAN THE RADIUS OF THE 
MARTIAN ATMOSPHERE 

IF (APOAP .LE. 3643.5) GO TO 80 
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45 WRITE ( 7, 120 ) I, PERAPE, APOAPE, SMAJE, PERDHR, DRAGE, TIME 

IF( APOAP .LE. APOAPF ) GO TO 60 
TIMTTL * TIMTTL*PERIOD 
50 CONTINUE 

IF (APOAP .GT. APOAPF) GO TO 85 

DETERMINE THE TIME TO TRAVEL THE HALF ORBIT FROM THE AEROBRAKING 
PERIAPSIS TO THE APOAPSIS. 

60 TIMEPA * . 5*PERI0D 

A DELTA-V BURN WILL BE PERFORMED AT THE APOAPSIS TO RAISE THE 
PERIAPSIS TO 500 km. DETERMINE THE PERIOD OF THE NEW ORBIT, AND 
THE TIME TO TRAVEL FROM THE APOAPSIS TO THE PERIARSIS. 

SMAJAP * . 5*( (500. *3393. 5) ♦ APOAP) 

PERDAP = 2. »PI*SGRT(SNAJAP**3/HU) 

TIMEAP » . 5*PERDAP 
IF (APOAP .EQ. APOAPF) THEN 
PERDF * PERDAP 
ELSE 

GO TO 70 
END IF 
GO TO 75 

BECAUSE THE FINAL APOAPSIS FROM AEROBRAKING IS LESS THAN THE DESIRED 
APOAPSIS, A DELTA-V BURN WILL HAVE TO BE APPLIED AT THE PERIAPSIS 
TO RAISE THE APOAPSIS SO THAT THE FINAL CIRCULAR ORBIT IS OBTAINED 
THE PERIOD OF THIS ORBIT, AND THE TIME TO COMPLETE ONE ORBIT (THUS 
FINALIZING THE CIRCULARIZATION OF THE ORBIT ABOUT MARS) IS DETERMINED 
70 SMAJF * 500. *3393. 5 

PERDF * 2. *PI*SQRT(SMAJF**3/MU) 

75 TIMTL * TIMTTL*TIMEPA*TIMEAP*PERDF 
TIMHR * ( TIMTL+TIME1 *TIME2 ) /3600. 


WRITE (7, •) 'TIME BREAKDOWN (hrs): ' 

WRITE(7, *) ' TIME TO INITIALIZE ORBIT: ', TIME1/3600. 

WRITE (7, • >' TIME TO TRAVEL FROM APOAP TO PERIAPs TIME2/3600. 

WRITE (7, *) ' TIME FOR AEROBRAKING PASSAGE: ' , TIHTTL/3600. 

WRITE (7, • ) ' TIME TO TRAVEL FROM PERIAP TO APOAPs ' , TIMEPA/3600. 

WRITE (7, •>' TIME TO TRAVEL FROM APOAP TO PERIAP: ' , TIMEAP/3600. 

WRITE (7, *) ' TIME FOR 1 ORBIT AFTER CIRCULARIZE: \ PERDF/3600. 

WRITE (7, * ) ' 

WRITE (7, • ) ' TOTAL TIME FOR AEROBRAKING PROCESS: TIMHR 

PRINT*, I 

PRINT*, 'AEROBRAKING TIME* ', (TIMTL«-TIMEl«-TIME2)/3600. 

PRINT*, 'PER TO APO* ', TIMEPA/3600. 

PRINT*, 'APO TO PER* ', TIMEAP/3600. 

PRINT*, 'ORBIT AFTER CIRCULARIZATION* ', PERDF/3600. 

GO TO 90 

80 WRITE(7, *)' AEROBRAKING IS NOT POSSIBLE FOR THIS PERIAPSIS' 

GO TO 90 

8S WRITE (7,*) 'FINAL APOAPSIS HAS NOT BEEL REACHED' 

130 FORMAT (3X, 'PERIAPSIS (ft) ', 6X, 'APOAPSIS (ft)',6x, 

1 ' SEMI -MAJOR AXIS', 5X, 'PERIOD (hrs)') 

140 FORMAT ( 2X, F15. 5, 5X, FIS. 5, 5X, FIS. S, 5X, Fll. 5) 

100 FORMAT (4X, 'PASS', 5X, 'PERIAPSIS', 7X, 'APOAPSIS', 7X, 

1 ' SEMI -MAJOR ' , 4X, 'PERIOD', 7X, 'DRAG',8X, 'PASSAGE') 

110 F0RMAT(3X, 'NUMBER', 7X, '(ft)', llx, '(ft)', llx, 'AXIS (ft)', 

1 5x, ' (hrs) ', 7x, ' (lb) ', 7x, 'TIME (min)') 

120 FORMAT (4X, 13, 5X, F13. 3, 3X, F13. 3, 3X, F13. 3, 3X, F7. 3, 3X, F10. 3, 4X, F8. 3) 

90 CLOSE (UNIT * 7) 

STOP 

END 


■ SUBROUTINE PARAMS(RP, RA, A, PERD) 

C THIS SUBROUTINE CALCULATES THE APOAPSIS AND PERIOD OF THE 

■ ELLIPTIC ORBIT, USING THE VALUES OF THE PERIAPSIS AND SEMI- 

| MAJOR AXIS FROM THE MAIN PROGRAM 

f THE PERIAPSIS, APOAPSIS, AND SEMI -MAJOR AXIS ARE IN km 
| THE PERIOD IS IN sec, THE GRAVITATIONAL PARAMETER IS <km*3/sec*2> 

REAL MU 

1 RA * 2. *A-RP 

■ PI * 3.141592654 

MU * 42656.0 

■ PERD * 2. *PI*SQRT(A**3/MU) 

I RETURN 

END 

| SUBROUTINE NTRSEC ( A, B, E, X, V ) 

C THIS SUBROUTINE CALCULATES THE POINTS OF INTERSECTION OF THE 

^ SPACE VEHICLE'S ELLIPTIC ORBIT AND THE ATMOSPHERE'S CIRCULAR 

K ORBIT, USING THE SEMI -MAJOR AXIS AND THE ECCENTRICITY FROM THE 

^ MAIN PROGRAM, AND THE SEMI -MINOR AXIS FROM SUBROUTINE PARAMS 



RADIUS - 250. >3393. 5 

XI « 2. *A»E 

X2A * 4. *A»»2»E»*2 

X2B * 4. *(B**2-RADIUS**2>A*»2*E»*2) *(i. -B»*2/A»*2) 

X2 > SORT ( X2A-X2B ) 

X3 ■ 2»<1. -B**2/A**2' 

THE INTERSECTION POINTS OF THE ELLIPTIC ORBIT ARE X AND Y 
X ■ (Xl-X2)/X3 

Y * SQRT(RADIUS**2-(X-A*E> **2) 

RETURN 

END 

SUBROUTINE SEGHNT(X, Y, A, E, SEG, PHI, AREA) 

THIS SUBROUTINE CALCULATES THE LENGTH OF THE SEGMENT (km) OF THE 
SPACE VEHICLE'S ELLIPTIC ORBIT BOUNDED BY THE MARTIAN ATMOSPHERE 
USING THE INTERSECTION POINTS, SEMI -MAJOR AXIS, AND ECCENTRICITY 
FROM THE MAIN PROGRAM 


C * 2. *Y 
D » X-A*E 

R * SQRT(D**2>Y<m»2) 

PHI IS THE ANGLE OF THE BOUNDED SEGMENT, AND AREA IS THE AREA 
OF THE BOUNDED PORTION OF THE ORBIT 
PHI * 2. *ATAN(C/(2. »D) ) 

SEG * R*PHI 
AREA > . 5*R*SEG 
RETURN 
END 



AEROBRAKE ANALYSIS 


ORIGINAL PAGE IS 
0F POOR Quaut? 


HALF- ANGLE FOR CONICAL AEROBRAKE (d*g) : 70.00000000 

DIAMETER OF THE AEROBRAKE (ft): 95.00000000 

SURFACE AREA OF THE AEROBRAKE (ft A 2): 7543.12402000 

HASS OF SPACE VEHICLE (lb): 200000.00000000 


ATMOSPHERIC CONDITIONS: 

PERIAPSIS FROM CENTER OF MARS (ft): 1. 14484910E*07 

DENSITY ( slug/ft A 3) : 2. 42937948E-10 

INITIAL ORBITAL PARAMETERS: 

PERIAPSIS (ft) APOAPSIS (ft) SEMI-MAJOR AXIS PERIOD (hrs) 

11448491.00000 120592192.00000 66020340.00000 24.12273 

APPROX TIME (hr) PRIOR TO AEROBRAKING? 48.79166030 


AEROBRAKE PROCEDURE: 


PASS 

PERIAPSIS 

APOAPSIS 

SEMI -MAJOR 

PERIOD 

DRAG 

NUMBER 

(ft) 

(ft) 

AXIS (ft) 

(hrs) 

(lb) 

1 

11448491.000 

113458000. 000 

62453248. 000 

22. 194 

388. 908 

2 

11448491.000 

107067800. 000 

59258144. 000 

20. 513 

386. 799 

3 

11448491.000 

101310824. 000 

56379656. 000 

19. 037 

384. 695 

4 

11448491.000 

96097224. 000 

53772856. 000 

17. 732 

382. 594 

5 

11448491.000 

91353344. 000 

51400916. 000 

16. 572 

380. 499 

6 

11448491.000 

87018344. 000 

49233416. 000 

IS. S3S 

378. 407 

7 

11448491.000 

83041360. 000 

47244924. 000 

14. 603 

376. 319 

8 

11448491. 000 

79379664. 000 

45414076. 000 

13. 762 

374. 235 

9 

11448491.000 

75997080. 000 

43722784. 000 

13. 001 

372. 155 

10 

11448491.000 

72862736. 000 

42155612. 000 

12. 30 8 

370. 079 

11 

11448491.000 

69950144. 000 

40699320. 000 

11.676 

368. 006 

12 

11448491.000 

67236472. 000 

39342480. 000 

11.097 

365. 937 

13 

11448491.000 

64701880. 000 

38075184. 000 

10. 565 

363. 872 

14 

11448491.000 

62329112. 000 

36888804. 000 

10. 075 

361. 810 

15 

11448491.000 

60103032. 000 

35775760. 000 

9.623 

359. 751 

16 

11448491.000 

58010348. 000 

34729420. 000 

9.204 

357. 695 

17 

11448491.000 

56039316. 000 

33743904. 000 

8. 815 

355. 642 

18 

11448491.000 

54179508. 000 

32813998. 000 

8. 453 

353. 592 

19 

11448491.000 

52421676. 000 

31935084. 000 

8. 115 

351. 545 

20 

11448491. 000 

50757564. 000 

31103028. 000 

7. 800 

349. 501 

21 

11448491.000 

49179764. 000 

30314128. 000 

7.505 

347. 459 

22 

11448491.000 

47681636. 000 

29565064. 000 

7. 229 

345. 419 

23 

11448491.000 

46257200. 000 

28852846. 000 

6.969 

343. 382 

24 

11448491.000 

44901056. 000 

28174772. 000 

6. 725 

341. 347 

25 

11448491.000 

43608300. 000 

27528396. 000 

6. 495 

339. 314 

26 

11448491.000 

42374508. 000 

26911498. 000 

6. 278 

337. 282 

27 

11448491.000 

41195628. 000 

26322060. 000 

6. 073 

335. 252 

28 

11448491.000 

40067988. 000 

25758238. 000 

5.879 

333. 224 

2q 

11448491. 000 

38988220. 000 

252183S4. 000 

5. 695 

331. 197 


PASSAGE 
TIME (min) 

11.691 
11.762 
11.835 
11.909 
11.983 
12.060 
12. 137 
12. 216 
12. 296 
12. 378 
12. 461 
12. 546 
12. 633 
12. 721 
12.811 
12. 902 
12. 996 
13.091 
13. 188 
13. 288 
13. 389 
13. 493 
13. 599 
13. 708 
13.819 
13. 933 
14. 050 
14. 169 
14. 292 


1 31 

11448491. 000 

36960228. 000 

S 32 

11448491.000 

36006576. 000 

1 33 

11448491.000 

35089896. 000 

34 

11448491.000 

34207976. 000 

■ 35 

11448491.000 

33358782. 000 

36 

11448491.000 

32540420. 000 

I 37 

11448491.000 

31751138.000 

1 38 

11448491.000 

30989304. 000 

39 

11448491.000 

30253406. 000 

1 40 

11448491.000 

29542030. 000 

I 41 

11448491.000 

28853862. 000 

42 

11448491.000 

28187662. 000 

1 43 

11448491.000 

27542282. 000 

44 

11448491.000 

26916638. 000 

1 45 

11448491.000 

26309712. 000 

. 46 

11448491.000 

25720544. 000 

47 

11448491. 000 

25148234. 000 

1 48 

11448491. 000 

24591922. 000 

49 

11448491. 000 

24050806. 000 

1 50 

11448491.000 

23524116. 000 

1 51 

11448491. 000 

23011118.000 

52 

11448491.000 

22511124.000 

| 53 

11448491. 000 

22023460. 000 

54 

11448491.000 

21547492. 000 

55 

11448491.000 

21082604. 000 

56 

11448491.000 

20628198. 000 

57 

11448491.000 

20183694. 000 

58 

11448491.000 

19748528. 000 

59 

11448491.000 

19322144. 000 

60 

11448491. 000 

18903988. 000 

61 

11448491. 000 

18493504. 000 

62 

11448491. 000 

18090140. 000 

63 

11448491.000 

17693322. 000 

64 

11448491. 000 

17302464. 000 

65 

11448491. 000 

16916940. 000 

66 

11448491.000 

16536060. 000 

67 

11448491.000 

16159161.000 

68 

11448491.000 

15785364. 000 

69 

11448491.000 

15413744. 000 

70 

11448491. 000 

15043186. 000 

71 

11448491. 000 

14672308. 000 

72 

11448491.000 

14299343. 000 

73 

11448491. 000 

13921898. 000 

74 

11448491.000 

13536531. 000 

75 

11448491.000 

13137854. 000 

76 

11448491.000 

12716375. 000 

AEROBRAKING 

COMPLETE 



TIME BREAKDOWN (hrs): 


TIME TO INITIALIZE ORBITs 

TIME TO TRAVEL FROM APOAP TO PERIAP: 

TIME FOR AEROBRAKING PASSAGE: 

TIME TO TRAVEL FROM PERIAP TO APOAP: 
TIME TO TRAVEL FROM APOAP TO PERIAP: 
TIME FOR 1 ORBIT AFTER CIRCULARIZE: 


24204360. 000 

5.355 

327. 146 

14. 546 

23727534. 000 

5. 197 

325. 122 

14. 679 

23269194. 000 

5. 048 

323. 098 

14.815 

22828234. 000 

4.905 

321. 075 

14. 955 

22403636. 000 

4.769 

319. 051 

15. 099 

21994456. 000 

4.639 

317. 027 

15. 247 

21599814. 000 

4.514 

315. 003 

15. 400 

21218898. 000 

4.395 

312. 978 

15. 557 

20850948. 000 

4. 282 

310. 953 

15. 720 

20495260. 000 

4. 172 

308. 925 

15. 887 

20151176. 000 

4.068 

306. 897 

16.061 

19818076. 000 

3.967 

304. 866 

16. 240 

19495386. 000 

3.871 

302. 833 

16. 426 

19182564. 000 

3. 778 

300. 797 

16.619 

18879102. 000 

3.689 

298. 758 

16.818 

18584518. 000 

3.603 

296. 715 

17. 026 

18298362. 000 

3. 520 

294. 669 

17. 242 

18020206. 000 

3. 440 

292. 618 

17. 467 

17749648. 000 

3.363 

290. 561 

17.701 

17486302. 000 

3. 288 

288. 500 

17. 946 

17229804. 000 

3.216 

286. 431 

18. 202 

16979808. 000 

3.146 

284. 356 

18. 470 

16735976. 000 

3. 079 

282. 273 

18.752 

16497992. 000 

3.013 

280. 182 

19. 048 

16265547. 000 

2. 950 

278. 081 

19. 360 

16038344. 000 

2.888 

275. 969 

19. 690 

15816093. 000 

2.829 

273. 846 

20. 039 

15598510. 000 

2. 770 

271.711 

20. 409 

15385317. 000 

2.714 

269. 561 

20. 804 

15176239. 000 

2. 659 

267. 395 

21.225 

14970998. 000 

2.605 

265. 213 

21.677 

14769316. 000 

2.552 

263.011 

22. 164 

14570907. 000 

2.501 

260. 788 

22. 690 

14375477. 000 

2. 451 

258. 540 

23. 261 

14182715. 000 

2. 402 

256. 266 

23. 886 

13992265. 000 

2. 354 

253. 961 

24. 574 

13803826. 000 

2.306 

251.622 

25. 337 

13616927. 000 

2.260 

249. 244 

26. 190 

13431118. 000 

2.213 

246. 820 

27. 155 

13245838. 000 

2. 168 

244. 344 

28. 260 

13060399. 000 

2.122 

241.805 

29. 547 

12873917. 000 

2. 077 

239. 192 

31.074 

12685194. 000 

2. 032 

236. 489 

32. 934 

12492511.000 

1.986 

233. 672 

35. 280 

12293172. 000 

1.938 

230. 708 

38. 391 

12082433. 000 

1.889 

227. 544 

42. 842 


36. 73029330 
12. 06136420 
472. 14114400 
0. 94430298 
1.02305233 
2. 05304074 


TOTAL TIME FOR AEROBRAKING PROCESS: 524.95318600 



APPENDIX B 

Optimization of Propallant Haas of a 
Solld-Propallant Rockat 
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★ ★ ** *★ kit if it kkkk ****** ****** kit kkkk kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 

MIKE LISANO 

AUEURN UNIVERSITY AEROSPACE ENGINEERING SENIOR DESIGN PROJECT 
UNPANNED MARS MISSION 


M ST AGES" 

PROGRAM TO OPTIMIZE THE PROPELLANT MASS OF A SOLID-PROPELLANT 
ROCKET LAUNCHING IN AN ARBITRARY GRAVITY FIELD FOR A GIVEN 
DESIRED ORBIT (PITCH PROGRAM INCLUDED) 
**************************************************************** 


C0MM0N/R0CK/TCCM9/DF/DIA/BR/GAMMA/RGAS/GC 
COMM ON/PLAN/GR A V/R/ A/OMEGA/ R AD/ ALAT 

C0MM0N/LAUN/V(610)/G(610)/H(61C)/AM(610)/DELV(610)/ACC(610) 

UNIT CONVERSION FACTOR GC (L3M-FT/LBF-SEC**2) 

GC*3 2.174 


* * ** ★* **** **************** ** ************ *************** **** ***** 


SPECIFICATIONS OF ROCKET: 

(COMBUSTION TEMPERATURE OF PROPELLANT/ DEGREES RANKINE) 
TCCM 8*6700. 

(DENSITY OF FUEL/ L8M/IN**3) 

DF*. 065 

(PROPELLANT CROSS-SECTION DIAMETER/ FT) 

FDIA*2. 91667 
DIA*FDIA*1<!.0 

(PROPELLANT BURN RATE/ IN/SEC) 

BR*0 .5 5 

(SPECIFIC HEAT RATIO OF BURNING PROPELLANT ( D£ FAULT : AI R ) ) 
GAMM A*1 • 4 

(PERFECT GAS CONSTANT OF BURNING PROPELLANT (DEFAULT : AIR) / 
FT-LBF/LBM-R) 

RG AS *5 3. 3 

(DEADWEIGHT RATIO OF BOOSTER) 

DWRA T* .1 2 

(MASS OF PAYLOAD TO BE CARRIED INTO 0R8IT/ L8M) 

AM F* 1000 .0 


**************************************************************** 

DATA FOR LAUNCH: 

(RAOIUS OF PLANET/ N MI) 

RAD* 1841 .05 

(LATITUDE OF LAUNCH SITE/ DEGREES) 

DL AT *0 .0 

ALAT*DLAT/53. 295773 

(ANGULAR VELOCITY OF PLANET/ RAD/SEC) 

OMEGA*. 000U7 

(DESIRED ALTITUDE/ N MI) 

ALT* 27C.0 
H F*ALT *608U. 
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(TOTAL HEIGHT/ FT) 

R*(RAD*608U.)+HF 

(SEMI-MAJOR AXIS OF DESIRED ORBIT/ N MI (DEFAULTS TO A*R) ) 

SMA*RAO+ALT 

A*$M A* 6080. 

(GRAVITATIONAL ACCELERATION ON THE SURFACE OF THE PLANET/ 

FT/ SEC**2) 

GSLR F*12.3J2 

(GRAVITATIONAL PARAMETER OF PLANET IN FT**3/SEC**2) 

GRAV*1 .5065E15 

**** ****** **************** ************** ************************ 


WR IT E( 6/ 1 1 U) 

,10 FORMAT (////3X/ 1 OPTIMIZATION OUTPUT:') 


SUBROUTINE "ROCKET" CALCULATES THE EXHAUST VELOCITY (C) AND TIME 
RATE OF CHANGE OF MASS (DMDT) OF THE ROCKET: 

CALL ROCKET(CzDMDT) 

WRIT E( 6/ 111) C/ DMDT 

11 FORMAT (///iX/* EXHAUST VELOCITY CF ROCKET *' / FI 2. 3/1 X/ ' FT/SEC ' / 
*///3Xz 'PROPELLANT CONSUMPTION RATE * ' / FI 2 .6/ 1 X/ • LBM/ SEC ' ) 


SUBROUTINE "SPEEDS" CALCULATES THE SPEED OF THE ROCKET DUE 
TO PLANETARY SPIN BEFORE TAKEOFF/ AND THE SPEED OF THE 
ROCKET IN THE DESIRED ORBIT: 

CALL SPEEDS(VSURFzVORB) 

WRITE(6/112)VSURF/ALT/V0RB 

12 F0RMAT(///3X/'VEL0CITY OF PLANET SURFACE * ' / FI 2. 3/ 1 X/ ' FT/SE C ' / 
*///3Xz 'DESIRED ORBIT ALTITUDE * ' / FI 2 . 3/1 X/ ' N MI'////3Xz 
★'VELOCITY OF ROCKET IN DESIRED ORBIT «' / FI 2. 3/1 X/ • FT/SEC ' ) 


SUBROUTINE "LAUNCH" CALCULATES THE TOTAL CHANGE IN SPEED AND 
CHANGE IN ALTITUDE OF A SINGLE STAGE ROCKET BEING LAUNCHED IN 
THE GRAVITY FIELD OF A GIVEN PLANET. AERODYNAMIC FORCES ON THE 
ROCKET HAVE BEEN NEGLECTED. THE INITIAL MASS OF THE ROCKET IS 
INCREMENTED FROM THE MASS REQUIRED FOR A TEN MINUTE BURN TO THAT 
REQUIRED TU REACH THE DESIRED ALTITUDE (HF) AND VELOCITY (V0R8). 

CALL LAUNCH(DWRAT/VSURF/G$URF/AMF/DMDT/C/HF/RAD/VORB) 

STOP 

END 

**************************************************************** 
SUBROUTINE ROC KET ( C/ DMDT ) 
COMMCN/ROCK/TCOMB/DF/DIA/BR/GAMMA/RGAS/GC 
- A*<3.14159j7*CDIA**2))/4. 

DMDT *A *D F*lsR 

C*SQRT ((2. *G AM MA*RGAS/ (GAMMA-1. ))*GC*TCOMB) 


S3 



I 

I 

I 

I 


RETURN 

END 

************************************ **************************** 
SUBROUTINE SPEEDS (VSURF/VORB) 

COPM ON/PLAN/ GRAV/R/ A/ OMEGA /R AD/ AL AT 
VStR F* CM EG A* RAD* 6080 •* COS ( AL AT) 

V0RB*SGRT(GkAV*((2./R)-(1./A))) 

RETURN 

END 

**************************************************************** 
SUBROUTINE LAUNCH (DWRAT/VSURF/GSUR F/AMF/DMDT/C/HF/ RAD/ VORB) 
COPMON/LAUN/V(610)/G(610)/H(61C)/AM(610)/D£LV(610)/ACC(610) 


JlO 

r 


TARGET BURN TIME (SEC). . . DETERMINES INITIAL MASS AM(1) 

TBURN*350.U 

TIRE INCREMENTS (SEC) 

D T=1 .0 

T*C. 0 

AM (1 )*((TBURN*DMDT)/(1 .O-DLR AT) ) + AMF 
AMFF*(AM(1; 'AMF)*DWRAT+AMF 
V(1) *VSURF 
DELV (1 )*0.U 
G ( 1 ) *GSURF 
H ( 1 ) *0.0 
DTHET*0. 0 
1*1 
1 * 1+1 
T*T+DT 

THET*DTHET/57. 29578 
AM (I)*AM(I-1)-(DMDT*DT) 

IF (AM(I) .LT.AMFF)GOTO 1900 

DELV(I)*C*L0G(AM(I“1)/AM (I ) ) - ( G ( 1-1 ) *C OS (THET) * DT ) 
V(I)*V(I-1)+DELV(I) 

H( I)*H(I-1 )+(< (V(I)+V(I-1) )/2.)*DT)*SIN(THET) 

6(I) a GSUHF*((RAD*6Q80.)**2)/(((RAD*6080.)+H(I))**2) 

ACC(I)*DfcLV(I)/(DT*32.174) 

PITCH PROGRAM (ARBITRARY/ FIVE STEP/ INITIAL THETA * 0 DEG./ 
FINAL THETA * 90 DEG. ) 

IF(H (D.GE.HF)GOTO 1900 
IF (H (I).GT.(HF*.4))G0T0 1800 
IF (H (I ).GT .(HF*.1 5))G0T0 1700 
IF(H (I).GT.(HF*.05))G0T0 1600 
IF(H(I).GT.(HF*.01))G0T0 1500 


«00 
i 60 0 

|oo 
1 800 


DTHET*A5.0 
GOTO ICO 
DTHET*65 .0 
GOTO ICO 
DTHE T* 80. 0 
GOTO 100 
DTHE T* 85 . 0 
GOTO ICO 



VEL«V(I-1) 

HT*H(I-1)/$080. 

AC*ACC(I“1 ) 

AMP* AM <1-1 ) 

PR*AMF /AM<1) 

TPR*AM(I-1 J/AM<1) 

EX C* AM (1-1 J-AMF 

WRITE<6/212)TBURN/AM(1)/T/VEL/HT/AC/AMM/PR/TPR/EXC 
FORMAT (/// JX/ ' DATA FOR SINGLE-STAGE ROCKET LAUNCH:*/ 

*///3 X/ *THE TARGET TIME OF SURN IS * / F8.3/1 X/ ' SEC ' / 

*/ / /3 X/ *THE INITIAL MASS OF THE ROCKET IS ' / FI 2 . 3/ 1 X/ * LBM ' / 

*///3 X/ *THE TOTAL TIME OF BURN IS * / F8.3/1 X/ ' SEC * / 

*///3 X/ *THE FINAL VELOCITY IS * >F1 2.3/1 X/ ' FT/ SEC * / 

*///3X/*THE FINAL ALTITUDE IS */ FI 2.4/1 X/ * N MI'/ 

*// /3 X/ *THE FINAL ACCELERATION IS ‘/F6.2/1X/ * GS'/ 

*/ / /3 X/ *THE FINAL MASS OF THE ROCKET IS * / FI 2 . 3/1 X/ ' LBM * / 

*/ / /3 X/ *THE DESIGN PAYLOAD RATIO OF THE ROCKET IS '/F9.6/ 

*///3Xz ‘ACTUAL FINAL MASS/IMTIAL MASS IS '/F9.6/ 

*// /3 X/ *THE EXCESS MASS AFTER FIRING IS * / FI 2 . 3/ 1 X/ * LBM * ) 

RETURN 

END 

************************************* *************************** 
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OPTIMIZATION OUTPUT: 

EXHAUST VELOCITY OF ROCKET * 8968.144 FT/SEC 

PROPELLANT CONSUMPTION RATE » 34.395610 LBM/SEC 

VELOCITY OF PLANET SURFACE * 783.551 FT/SEC 

DESIRED ORBIT ALTITUDE « 270. OCO N MI 

VELOCITY CF ROCKET IN DESIRED ORBIT * 10833.867 FT/SEC 

DATA FOR SINGLE-STAGE ROCKET LAUNCH: 

THE TARGET T IhE OF BURN IS 250. OCO SEC 

THE INITIAL MASS OF THE ROCKET IS 10771.480 LBM 

THE TCTAL TIME OF BURN IS 250. 000 SEC 

THE FINAL VELOCITY IS 14047.491 FT/SEC 

THE FINAL ALTITUDE IS 210.3703 N MI 

THE FINAL ACCELERATION IS 4.35 GS 

THE FINAL MASS OF THE ROCKET IS 2172.578 LBM 

THE DESIGN PAYLOAD RATIO OF THE ROCKET IS 0.092838 

ACTUAL FINAL MASS/INITIAL MASS IS 0.201697 

THE EXCESS MASS AFTER FIRING IS 1172.578 LBM 
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OPTIMIZATION OUTPUT 


EXHAUST VELOCITY OF ROCKET * 8968.144 FT/SEC 

PROPELLANT CONSUMPTION RATE * 34.395610 LBM/SEC 

VELOCITY OF PLANET SURFACE * 783.551 FT/SEC 

DESIRED ORBIT ALTITUDE » 270. OGO N MI 

VELOCITY OF ROCKET IN DESIRED ORBIT * 10833.867 FT/SEC 

DATA FOR SINGLE-STAGE ROCKET LAUNCH: 

THE TARGET TIME OF BURN IS 350.000 SEC 

THE INITIAL MASS OF THE ROCKET IS 14680.072 LBM 

THE TOTAL TIME OF BURN IS 336.000 SEC 

THE FINAL VELOCITY IS 13328.603 FT/SEC 

THE FINAL ALTITUDE IS 268. 7500 N MI 

THE FINAL ACCELERATION IS 3. 03 GS 

THE FINAL MASS OF THE ROCKET IS 3123.147 LBM 

THE DESIGN PAYLOAD RATIO OF THE ROCKET IS 0.068120 

ACTUAL FINAL MASS/INITIAL MASS IS 0.212747 

THE EXCESS MASS AFTER FIRING IS 2123.147 LBM 
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OPTIMIZATION OUTPUT 


EXHAUST VELOCITY OF ROCKET = 8963.144 FT/SEC 

PROPELLANT CONSUMPTION RATE * 34.395610 LBM/SEC 

VELOCITY CF PLANET SURFACE * 733.551 FT/SEC 

DESIRED ORBIT ALTITUDE * 270. CCO N MI 

VELOCITY CF ROCKET IN DESIRED CR3IT a 10833.867 FT/SEC 

DATA FOR SINGLE-STAGE ROCKET LAUNCH: 

THE TARGET TIrE OF BURN IS 450. OCO SEC 

THE INITIAL MASS OF THE ROCKET IS 18588.664 L3M 

THE TCTAL TIME OF BURN IS 394.000 SEC 

THE FINAL VELOCITY IS 10955.684 FT/SEC 

THE FINAL ALTITUDE IS 268.4356 N MI 

THE FINAL ACCELERATION IS 1.37 GS 

THE FINAL MASS OF THE ROCKET IS 5036.794 LBM 

THE DESIGN PAYLOAD RATIO OF THE ROCKET IS 0.053796 

ACTUAL FINAL TASS/INITIAL MASS IS 0.270961 

THE EXCESS MASS AFTER FIRING IS 4036.794 LBM 


SS 
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ALT 




FINAL VELOCITY, FT /SEC 


ORIGINAL PAGE IS 
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FIGURE B. 2 

FINAL VELOCITY VS. "DESIGN BURN TIME" 



INITIAL PROPELLANT MASS/BURN RATE, SEC 

■< 10»*+Q2 
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PAYLOAD RATIO 


ORIGINAL PAGE IS 
0F POOR QUALITY 


FIGURE B. 3 



0.000 1.000 2.000 3.000 - 1.000 5.000 6.000 7.000 6.000 


INITIAL PROPELLANT/BURN RATE, SEC 

■< 10»*+Q2 
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RATIO 



0. 000 


I I I I 

0 . 000 1 . 000 1 

INITIA 







FINAL ACCELERATION, G* S 


PAGE ls 
Of POOR OUAUT? 


FIGURE B. 5 


FINAL ACCELERATION VS. "DESIGN BURN TIME" 



INITIAL PROPELLANT MASS/BURN RATE, SEC 

< i 0-**+02 
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APPENDIX C 

Calculation of Stagnation Taaparatura on Mara Lander 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


AppIy .tht aathad g< Mlaalal 4t.fr l» 


Fro* Sttlan’fl Lav f 


T„ 


i **» » • 
6’ffn 


.ts 

) 


vhara T« la tht ttaparatur* at th# stagnation point of tha body 
(dagraas Rankina) 

qm is tha radiativa haat flux fros air to body 
(Btu/ft* /sac) 


€ is tha asissivity of tha fluid 

is tha Stafan-Boltzaan constant 
<.481 E-12 Btu/ft* /sac/ °R) 


For tha atsosphara of Kars tha asissivity is assusad to ba 
approxlsataly tha sasa valua as tha asissivity of tha Earth's at 
aosphara, or 

€ - 0.8 

For tha flov-body systas to ba in aqulllbrius, tha haat 
radlatad to tha body sust aqual tha haat eonvactad fros tha body 

to tha flows 

• • 

q«*t * q« • v 

An aspirical forsula fros Nicolai allows tha calculation of 
tha convactiva haat flux: 

q*..* - lSM^j^tjggj'-cos A 

vhara o » fraastraas danslty (slugs/ft 1 > 

J 60 

u ■ fraastraas valocity (ft/sac) 

R * radius of curvatura of tha nosa (ft) 

A » svaap angla of wing loading adga (zero dagraas) 
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At an altitude of lOO nautical ailaa tha density of the 
Martian ataosphere is deterainedt 

_ * P >0 ° 

Rcoj’Tw. 

vhere R ■ 35.10 ft-lbf/(lbs-*R) 

Co x 

- 324.6 °R 

P iook .^ • 0.0042837 lbf /ft* 

The values of T and P were obtained froa Viking data. 

The velocity at 100 nautical ailes vas estiaated by aasuaing 
that the specific total aechanical energy of the lander at 100 
nautical ailes is the sane as the specific total aechanical 
energy at 270 nautical ailes. This assuaption iaplies that no 
vork is done by drag forces, vhlch vould decrease the total 
energy of the lander. Thus the velocity estiaate is high: 
u Mj » 11,744 ft/sec 

too 

This estiaate leads to a convective heat flux of 
qsaav • 2.7199 Btu/f t* /sec * q B • • 
and a stagnation taaparature of 
T. - 1630 °R 
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APPENDIX D 

Optimization of Racovary Syataa Haaa 
For Mara Landars 
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RECOVERY SYSTEM SUMMARY 


Weight Statement 


Parachute System Weight 

8 

224 
69 

76 

391 X 3 
59 

SRM Weight 


Propellant 2243 

Case Weight. 396 


Pilot Chute. 

Drogue Chute 

Misc. Straps 

Main Parachute 

Support Structure 

Main Parachute. ......... 

Fittings and Flotation.. 


1610 


+ 2639 


Recovery System Total Weight 


4249 


Basic Vehicle Weight ' + 58500 

Total Vehicle Weight (Reentry) 62749 


Performance Characteristics 


Parachute System SRM 


1 pilot @ 11.7 ft. diameter ISP 260 sec 

1 drogue @ 54.4 ft. diameter mass fraction 0.850 

3 mains @ 129.7 ft. diameter thrust 100000 lbf 

75 fps terminal velocity burn time 6.04 sec 

impact velocity 9.92 fps 

constant velocity 

falling height 5.1 ft 

total distance fallen. 210.2 ft 

maximum deceleration.. 1.6 g's 
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MASSES VERSUS TERMINAL VELOCITY 

5000 -| -B- SRM MASS (LB) 

v * — — 17 PARACHUTE MASS (LB) 

-9- TOTAL MASS (LB) 



iv 
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o 
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o 

c 


rO 

CM 



(an) S355V1A! 
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MAIN CHUTE TERMINAL VELOCITY (FT/SEC) 



FI0URE D - 2 MASSES VERSUS IMPACT VELOCITY 

5000 -i -e- SRM MASS (LB) 

" _____ PARACHUTE MASS (LB) 

-v- TOTAL MASS (LB) 



(an) s3SSviai 



o 

O 

c> 

o 

o 

o 

o 

a 

o 
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a 


ro 

CM 
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IMPACT VELOCITY (FT/SEC) 


CD 



(an) sassvw 
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CONST. VEL. FALLING HT. (FT.) 


